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基于 Perron 补 的 2 矩阵 最 小 特征 值 界 的 估计 


杨志明 
(甘肃 联合 大 学 数 信 学 院 ， 兰 州 730000) 
摘 要 : 本 文 给 出 了 估计 不 可 约 2 矩阵 的 最 小 特征 值 上 下 界 的 一 种 简单 方法 ， 即 以 矩阵 的 广义 Perron 补 
为 基础 ， 将 不 可 约 2Z- 和 矩阵 4 = 5 一 吾 的 最 小 特征 值 问题 化 为 广义 Perron 补 已 _。B)(4/4a) 的 
最 小 特征 值 问 题 ， 然 后 利用 矩阵 范 数 的 性 质 导 出 了 4 的 最 小 特征 值 界 的 估计 式 ， 同 时 也 给 出 了 
非 负 不 可 约 矩 阵 殖 的 谱 半径 的 一 种 简单 估计 式 . 
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1 引言 


2 矩阵 是 比 M- 和 矩阵 更 广泛 的 矩阵 . 工程 实践 中 的 许多 振动 问题 ， 如 桥梁 、 建 筑 物 的 振 
动 、 机 械 机 件 、 飞 机 机 副 的 振动 及 稳定 性 分 析 问 题 常 需要 对 其 最 小 特征 值 的 界 进行 估计 . 
文献 中 利用 相似 变换 得 到 了 一 种 同步 计算 最 小 特征 值 及 特征 向 量 的 方法 ， MeyerD 和 Lull 先 
后 提出 了 非 负 不 可 约 矩 阵 的 Perron 补 及 广义 Perron 补 的 概念 . 本 文 以 这 些 概念 为 基础 ， 通 
过 2G- 抢 阵 与 非 负 移 阵 及 其 广义 Perron 补 的 关系 来 研究 Z- 和 矩阵 的 最 小 特征 值 的 上 下 界 问题 . 

2- 和 矩阵 就 是 非 对 角 线 上 的 元 素 小 于 等 于 零 的 矩阵， 可 表 为 人 4= sT 一 已 ， 其 中 媚 之 0，T 为 
单位 矩阵 ， 而 且 若 4 为 不 可 约 矩阵 ， 则 B 为 非 负 不 可 约 和 矩阵 . 

若 分 别 用 w(A) 和 p( 吾 ) 表示 拢 阵 4 的 最 小 特征 值 和 非 负 不 可 约 扼 阵 吾 的 谱 半径 ， 则 


w(4) =s 一 p(B)， 上 
引 理 1 四 “ 设 4 = (asy)wxn 为 非 负 不 可 约 矩 阵 ， 且 m = > aij (1 < i<n)， 则 
minr;=7(A) < p(A) < R(A) = maxrs. (2) 
引 理 2 四 ”车 矩阵 A 的 范 数 | 4 < 1， 则 了 一 和 4 非 奇异 ， 且 


1 
I— 4) "|< . 
I ) | < 7a 
男 设 (n) = {1,2,… ,n}，aQ 为 (n) 的 非 空 有 序 子 集 ，B8 = (n) \ a， 则 有 : 
定义 1 四 设 4 = (ai) 是 谱 半径 为 p(4) 的 m 阶 非 负 不 可 约 矩 阵 ， 对 a C (n)， 用 A[a, 8] 
表示 元 素 为 uij 的 子 矩阵 , 其 中 ie a j € 8，4[a] 等 价 于 A[a,a]， 则 [a] 的 广义 Perron 补 定 
义 为 
P(A/Alal) = A[B]+ AlB,al(tI — Ala]) Alo, BJ, t>p(Alal). 
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当 t = p(4) 时 ， 称 之 为 Perron 补 ， 并 记 为 


P(A/Ala]) = AIB] + AlB, a](p(A)I — Ala]) ’ Ala, Al. 
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引 理 3l3] 设 A 是 谱 半 径 为 p(A) 的 非 负 不 可 约 和 矩阵 ， 则 对 任意 t > p(4[a)，4[a] 的 广 


义 Perron 补 P(A/A[al) 也 是 非 负 不 可 约 和 矩阵 ， 且 
(a) 当 t=p(4) 时 , p(R(A/Alal)) = p(P(A/Alal)) = p(4); 
(b) p(B(A/A[al])) 是 关于 t 的 严格 递减 函数 . 

2 w(PB(4/A[lal)) 与 p(B) 之 间 的 关系 
对 QC (n), B= (n) \a， 为 叙述 方便 ， 记 


Blal = Bo, Bla, A] = Bag, BlB,al = Bpo, BIB] = Bg, 
“ 名 
Bpa Bp 


[stI_B, _B 
A=sI_- B=|° a oh eB 
—Bpga sT— Bg 


且 不 妨 设 


则 


于 是 ，Ba 的 Perron 补 
P(B/Bo)= Bp+ Bpa(p(B)I- Ba) Bag, p(B)> p(Bo), 
而 4。 的 广义 Perron 补 
Pt(A/Aa) = sI -Ba+(-Bga)[(s —t)I—(sIT— Bo)] (~-Bos) 
= sT— [B+ BpaltT -~ Ba) ' Bog] = sT ~— P(B/Bo), 


其 中 s 一 +> p(Aa)， 由 (3) 式 即 得 : 
定理 1 设 4=sT-- 态 ( 百 >0) 为 不 可 约 G- 拖 阵 ， 且 s >t+p(4o)， 则 


w(Ps_1(A/Aa)) =s—p(P(B/Bo)). 


由 (1) 式 及 引 理 3(a)， 显 然 有 : 
推论 1 设 和 A= sI 一 B 为 不 可 约 2- 和 矩阵 ， 且 s > p(B) 十 p(Aa)， 则 


w (Ps_p(B)(A/Aa)) =s— p(B)=w(A). 
特别 地 ， 若 s = 2p(B)， 则 


w (PB)(A/Ao)) = w(4) = p(B). 


(3) 


这 就 是 说 ， 计 算 不 可 约 2 矩阵 4 的 最 小 特征 值 , 只 需 计算 广义 Perron 补 Psp)(4/ 4a) 的 
最 小 特征 值 ， 但 此 时 却 需 要 知道 p(B) 的 值 ， 因 而 实际 操作 中 是 行 不 通 的 。 事实 上 ， 我 们 可 以 


利用 推论 1 的 结论 来 估计 p(B) 的 界 ， 进 而 确定 w( 4) 的 界 . 
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3 不 可 约 2G- 和 矩阵 最 小 特征 值 的 界 
首先 ， 由 (1) 式 及 引 理 1 可 知 ， 不 可 约 Z- 抢 阵 4 的 最 小 特征 值 满足 
— R(B) < w(A) < s—"7(B). (4) 


另 由 推论 1 知 ，s p(B) 是 

Po(B)(4/4[a) = sI- [Be+ Bpa(p(B)I — Ba) Bog] 
的 一 个 特征 值 ， 所 以 存在 x 关 0， 使 

[sT— (Bg + Bgpa(p(B)I— Ba) Bop)]jz=(s-p(B))z 


成 立 ， 两 边 取 范 数 并 由 范 数 的 性 质 ||C + DD < Cl + | 刀 |， 其 中 C, 万 为 同型 矩阵 ， 可 得 


ls — p(B)| 1 
< |lsT- (Bs + Boa(p(B)T- Ba) Bog))| :zl 
< [llsIl+ |Bell+ lBgall: (p(B)IT — Ba) -Bael] -lzll, (5) 


车 p(B) > Ba。||， 则 由 引 理 2 知 矩 阵 p(B)I 一 B。 非 奇异 ， 且 


|(p(B)T- Bo | = zm |- BD) | i 


又 因为 s > p(B) 十 p(A。) > p(B), ||z|| > 0， 所 以 (5) 式 可 化 为 


-p(B) < s+ Bol + SE, 
即 
p(B) - (Bcl ~ Bal)p(B) + (IBaoll :Bagll -1Boll :|Bgl) > 0. 
由 此 得 


p(B) > 3[ (Bal ~ Bol) + VIBol + Bol)’ -4lBooll -1Boal | Eh (0 
完全 类 似 地 ， 由 范 数 的 性 质 咱 Cl - |PI < lc + 吃 上 | 又 可 得 到 

p(B) < 3[ (Bol + Bol) + V(IBol -Bol) +4IB6ol Bo | Eb 
综 上 可 得 以 下 定理 : 


定理 2 设 A= sI 一 B 为 不 可 约 2Z- 和 矩阵 ， 则 当 p(B) > B64 时 , 三 < p(B) < ts。. 进而 
有 s 一 t2 <w(A)<s-t. 
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4 ”数值 例子 与 界 的 改进 
例 1 对 2 矩阵 4= sT 一 B， 其 中 








DD 直人 口 
OD Oo 呈 0 WW 
一 器 WW OO 
一 
Om %% ~ 
OA 记 一 AA % 


显然 r( 互 ) = 6,，R(B) = 34， 故 由 (4) 式 可 得 s 一 34 < w(4) < s 一 6. 现 取 a = {6},， 则 p = 
{1,2,… ,5}， 即 将 矩阵 妃 分 块 为 





若 用 矩阵 的 行 范 数 上 ll。， 则 Ballw = 0, 1Bogllo = 6, Beallo = 8, 1Ball~ = 27， 从 而 
由 (6), (7) 两 式 得 


1 1 
一 1.9134 = 引 —27+ V27?—4x6x8|]<p(B)< 3[27+ V272 +4x6x8]= 28.6740, 


进而 有 s 一 28.6740 < w(4) < s 十 1.9134. 

显然 ， 这 里 得 到 的 p(B) 的 下 界 = 一 1.9134 没 有 (2) 式 给 出 的 下 界 r(B) = 6 精确 ， 事实 
上 ， 由 (6) 式 给 出 的 p(B) 的 下 界 在 许多 情况 下 是 无 法 使 用 的 ， 因 为 被 开 方 数 很 可 能 小 于 零 ， 
但 上 界 却 更 接近 于 p(B) (= 20.5044). 正 因为 这 样 ， 实 际 应 用 时 常用 (7) 式 估 计 p(B) 的 一 个 上 
界 ， 而 用 下 面 的 定理 确定 一 个 更 为 精确 的 下 界 . 


定理 3 设 
Bac Bp 
为 非 负 不 可 约 和 矩阵 ， 有 p(B) > Bll; t2 由 (7) 式 给 出 ， 则 p(B) > p(P,(B/B.o,)). 
证 明 由 (7) 式 知 p(B) < t。， 故 由 引 理 3， 得 
p(B) = p(PLB)(B/Bu)) > p(P(B/Bo)). 
推论 2 设 A=sI 一 B 为 不 可 约 2- 和 矩阵 ， 则 当 p(B) > 上 Bo|| 时 ， 有 


s—to <w(A)<s-p(P,(B/Bo)) =w(P,-1(B/Bo)). 
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如 例 1 中 取 a = {6} 时 ， 由 前 面 的 计算 知 t。 = 28.6740. 进而 有 p(P,(B/Ba)) = 20.0662， 
即 Z- 矩 阵 4 的 最 小 特征 值 满足 s 一 28.6740 < w(4) < s 一 20.0662. 这 显然 比 s 一 34 < w(4) < 
s 一 6 精确 得 多 . 

本 文 给 出 了 估计 不 可 约 GZ- 矩阵 4 = sIT -BB 的 最 小 特征 值 上 下 界 的 一 种 简单 方法 ， 其 实质 
还 在 于 估计 非 负 不 可 约 矩 阵 BB 的 谱 半 径 ， 即 


p(P,(B/Bo)) < p(B) < ts, (8) 


且 由 (8) 式 给 出 的 下 界 非 常 接近 于 p(B) 的 真 值 ， 如 例 1 中 非 负 不 可 约 矩 阵 妃 的 谱 半 径 p(B) = 
20.5044， 而 由 (8) 式 得 到 p(B) > 20.0662. 正 因 为 如 此 ， 本 文 的 方法 也 可 用 于 估计 非 负 不 可 
约 矩 阵 谱 半径 的 上 下 界 ， 尤 其 是 下 界 . 表 1 给 出 了 由 Matlab 随机 生成 的 元 素 介 于 0 ~ 9 之 间 
的 n 阶 非 负 不 可 约 矩 阵 吾 (=fix(10*rand(n)))， 当 n = 10, 30, 50 时 ， 由 (2), (8) 两 式 得 到 的 上 下 
界 ， 可 以 看 出 ， 本 文 的 方法 能 够 在 很 大 程度 上 改进 经 典 的 Frobenius 不 等 式 (2). 


表 1: 非 负 不 可 约 抵 阵 的 谱 半 径 估计 


公式 (2) 公式 (8) | 
p(B) 
r(B) R(B) p(Pa(B/Ba)) | 如 


10 39.0000 57.0000 46.6827 (a = {6}) 56.1127 47.6501 
30 92.0000 161.0000 134.2035 (a = {17}) 157.0000 135.1406 
50 183.0000 269.0000 226.0547 (a = {48}) 267.0671 226.8915 
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Estimating the Bounds of the Smallest Eigenvalue of Z-matrix with 
Perron Complement 


YANG Zhi-ming 


(College of Mathematics and Information, Gansu Lianhe University, Lanzhou 730000) 


Abstract: In this paper, we present a simple method to estimate the lower and upper bounds for the 
smallest eigenvalue of irreducible Z-matrices, the method is based on the generalized Perron comple- 
ment. For the smallest eigenvalue problem of the irreducible Z-matrix 4 = s7 — B, we convert it into 
the smallest eigenvalue problem of a generalized Perron complement. Then we utilize the properties of 
matrix norms and obtain the estimation of the bounds for the smallest eigenvalue of A. Moreover, we 
give a simple estimation for the spectral radius of a nonnegative irreducible matrix. 
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